Introduction

Animals come in many different shapes and sizes and they can be seen interacting with each other in many forms. One type of interaction between organisms is the predator-prey interaction which allows us to generate food webs. Cohen et al(1993) was able to determine that most predators that we see are bigger than the prey that they eat and there is a positive correlation for size increase in predators when prey size increases. Animals that eat at a higer trophic levels don’t neccessary have larger predator-prey ratios as discovered by Tucker and Rogers(2014) for marine and terrestrial mammals. There have been many other experiments that compare the mass of predators to the mass of prey and they overall agree that there is a relationship between predator-prey mass. However, are there other factors that can affect this relationship?

From the Marine predator and prey body sizes dataset created by Barnes et al (2008), we want to show that location and feeding behavior are just two of the aforementioned factors that could influence this interaction. The dataset contains a vast amount of information that comes from 27 marine locations with 93 types of predators and prey.

Fig.1: *Shows the log of Predator mass vs log of Prey Mass*

Fig.1: Shows the log of Predator mass vs log of Prey Mass

Fig.1: *Shows the log of Predator mass vs log of Prey Mass*

Fig.1: Shows the log of Predator mass vs log of Prey Mass

Fig.2: *Shows the log of Predator mass vs log of Prey Mass by Specific Habitat*

Fig.2: Shows the log of Predator mass vs log of Prey Mass by Specific Habitat

Fig.4

Fig.4

Methods

Data was obtained from the “Marine predator and prey body sizes”, created by Barnes et al(2008).; joint venture generated by 17 different international institutions. As such, we have data from global locations on organisms and their environment. We focused on data on environmental variables to see which ones suggest a correlation with predator mass and prey mass. Our variables of interest were selected to be: specific habitat location and feeding behaviour.

Importing the dataset into R, we used the package ggplot2 to visualize the patterns and trend in the aforementioned variables. Three graphs were generated from the data. The first figure compares predator mass and prey mass. The second figure is a map graph representing the locations of all the 27 international marine locations where our data was obtained. Using the facet feature of ggplot2, the third figure visualized the affect of feeding behaviour and the specific location on the relationship of predator-prey mass ratios.

To determine which variables are affecting the slope of predator-prey mass relations, we generated several models. (using glm; and lmr; and mn <- talk about these later). Every model assumes that the predator mass are the independent variable, and the prey mass are the dependent variable. The first model was a control model; assuming that there was nothing affecting the interaction of predator-prey mass. The second model incorporates the effect of specific habitat on the intercept and slope of the predator-prey mass ratio. Fourth model incorporates the effect of feeding behaviour on the intercept and slope of the predator-prey mass ratio. The third model incorporates the additive effect of both specific habitat and feeding behaviour. The sixth model incorporates the mixed effect of habitat on feeding behaviour in the predator-prey mass interaction. The seventh model incorporates the mixed effect of feeding behaviour in the mass interaction. The last model incorporates for “everything” (<- need more details).

Graphs

## We recommend that you use the dev version of ggplot2 with `ggplotly()`
## Install it with: `devtools::install_github('hadley/ggplot2')`

Results

The first graph shows a positive correlation between the log of predator mass and the log of prey mass; visualized via a linear regression line. As the dataset indicated, prey mass was determined off predator mass; ergo prey mass is the dependent variable.

Second graph is a map-graph representing the locations of the sample locations (MORE TO COME)

Third graph is the one that details how feeding behaviour and specific habitat affect the interaction of predator-prey mass. We can generalize the feeding behaviours of the predators into two categories: specialist feeders (the predators categorized as piscivorous, insectivorous, and herbivorous) and generalist feeders (categorized only has predacious). All the graphs under the specialist feeders showed a mild to strong positive correlation between predator and prey mass. In the piscivorous column, of the ten total plots, 8 of them show a positive correlation, 2 of them show neutral correlation In the planktivourous column, of the six plots that have a linear regression, 5 of them are positive in correlation, the last one shows a neutral correlation. On the other hand, within the graphs under generalist feeders, there were a variety of trends; some showed positive correlations, but also some interactions were neutral and even negative in correlation. *There are 12 plots in the column under predacious. 6 of the plots show an indication of positive trend, 5 show a neutral correlation, and one show a negative correlation.

Models

(Intercept) log(Predator_mass) class df logLik AIC delta weight
lm12 -5.869452 0.8529002 lmerMod 12 -72055.08 144134.2 0.000000 0.690088197
lm11 -5.375857 0.8361708 lmerMod 9 -72058.88 144135.8 1.613665 0.307965271
lm10 -5.910870 0.7056486 lmerMod 9 -72063.95 144145.9 11.741540 0.001946532
lm9 -6.500780 0.7208025 lmerMod 5 -72922.91 145855.8 1721.671423 0.000000000
null -7.027674 1.0048348 glm 3 -88088.81 176183.6 32049.463144 0.000000000

Table.1: Results of our models

## refitting model(s) with ML (instead of REML)
Df AIC BIC logLik deviance Chisq Chi Df Pr(>Chisq)
null 3 176183.6 176209.0 -88088.81 176177.6 NA NA NA
lm9 5 145850.5 145892.7 -72920.23 145840.5 30337.158841 2 0.0000000
lm10 9 144146.1 144222.2 -72064.07 144128.1 1712.310086 4 0.0000000
lm11 9 144136.5 144212.6 -72059.25 144118.5 9.654834 0 0.0000000
lm12 12 144135.6 144237.1 -72055.82 144111.6 6.854474 3 0.0766837

Discussion

Software

The software we used to create this document is the statistical Analysis program R. The packages used to create this data are:

  1. tidyverse
  2. dplyer
  3. readxl
  4. broom
  5. lme4
  6. MuMIn
  7. knitr
  8. plotly

References

  1. Cohen, J. E., Pimm, S. L., Yodzis, P., & Saldana, J. (1993). Body sizes of animal predators and animal prey in food webs. Journal of Animal Ecology, 62, 67-78. doi:10.2307/5483

  2. C. Barnes, D. M. Bethea, R. D. Brodeur, J. Spitz, V. Ridoux, C. Pusineri, B. C. Chase, M. E. Hunsicker, F. Juanes, A. Kellermann, J. Lancaster, F. M?nard, F.-X. Bard, P. Munk, J. K. Pinnegar, F. S. Scharf, R. A. Rountree, K. I. Stergiou, C. Sassa, A. Sabates, and S. Jennings. 2008. Predator and prey body sizes in marine food webs. Ecology 89:881.http://www.esapubs.org/archive/ecol/E089/051/metadata.htm

  3. Tucker MA, Rogers TL. 2014, Examining predator–prey body size, trophic level and body mass across marine and terrestrial mammals. Proc. R. Soc. B 281:20142103., http://dx.doi.org/10.1098/rspb.2014.2103